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[ We consider a discrete time random walk in one dimension. At each time step the walker jumps 

^vq . by a random distance, independent from step to step, drawn from an arbitrary symmetric density 

function. We show that the expected positive maximum E[M n ] of the walk up to n steps behaves 



3 



asymptotically for large n as, E[M„]/(T = yj / 2n/TT + j + C^n" 1 / 2 ), where a 2 is the variance of the 
step lengths. While the leading y/n behavior is universal and easy to derive, the leading correction 
term turns out to be a nontrivial constant 7. For the special case of uniform distribution over [— 1, 1], 
CofTmann et. al. recently computed 7 = —0.516068 ... by exactly enumerating a lengthy double 
series. Here we present a closed exact formula for 7 valid for arbitrary symmetric distributions. We 
also demonstrate how 7 appears in the thermodynamic limit as the leading behavior of the difference 
«4 ■ variable _E[M n ] — _E[|a; n |] where x n is the position of the walker after n steps. An application of 

^ ' these results to the equilibrium thermodynamics of a Rouse polymer chain is pointed out. We also 

generalize our results to Levy walks. 

PACS numbers: 02.50.-r, 89. 75. He, 89.20.Ff 
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Brownian motion is perhaps one of the most widely studied subjects in classical physics. A Brownian walker moves 
""O ■ in continuous space and time which leads to a great simplification: one can write down a simple differential equation, 
the famous diffusion equation, that governs the time development of the probability density of the position of the 
walker. Subsequently, many other involved properties of the Brownian motion, such as its first-passage probability 
through a given point, the distribution of the maximum displacement of the walker up to a given time etc. can be 
calculated analytically relatively easily [1]. In contrast, the related problem of a random walker that hops only at 
discrete time steps in a continuous space is not so straightforward, even in one dimension [2,3]. A classic example of 
such a walk can be found in bacterial chemotaxis, where a bacteria, in search of food, jumps from one position to 
' another at discrete time steps [4]. Another famous example of such a walk occurs in the Rouse model of a polymer 
t-H , chain that consists of monomers or beads connected by harmonic springs [5]. Many other examples can be found in 
Refs. [6-8]. While the asymptotic properties of this discrete hopper, after a sufficiently large number of steps and 
given that the variance in step sizes is finite, are correctly described by the continuous time diffusion equation [8], 
there are many interesting finite size effects that can not be captured by the diffusion equation. The difficulty that 
arises in dealing with a finite number of steps is due to the fact that the probability density of the random hopper 
usually satisfies an integral equation which is technically much harder to solve than a differential equation. In this 
paper, we study analytically one such finite size effect, namely the behavior of the expected maximum position of a 
discrete time hopper in one dimension. We will show that even this relatively simple problem has rather interesting 
finite size behavior. 

We consider a discrete time random walker moving on a continuous line. The position x n of the walker after n 



(J , steps evolves for n > 1 via, 

•En "b £n (1) 



' starting at xo = 0, where the step lengths £„'s are independent and identically distributed (i.i.d) random variables 
with zero mean and each drawn from the same probability distribution, Prob(£„ < x) = J** f(y) dy, f(x) being the 
normalized symmetric probability density. Let M n denote the positive maximum of the random walk up to n steps, 

M„ = max(0, ^i,ar 2 , . . . ,x n ). (2) 

We are interested in the asymptotic large n behavior of the expected maximum E(M n ). 

This question recently arose in the context of a packing problem in two dimensions where n rectangles of variable 
sizes are packed in a semi-infinite strip of width one [9,10]. It was shown in Ref. [10] that for the special case of the 
uniform jump distribution, f(x) = 1/2 for — 1 < x < 1 and f(x) = outside, for large n, 
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E[M n ) = J— - 0.297952 . . . + O^ 1 / 2 ). 

V 07T 



(3) 



The leading y/n behavior is easy to understand and can be derived from the corresponding behavior of a continuous 
time Brownian motion after a suitable rescaling [10]. However, the leading finite size correction term turns out to 
be a nontrivial constant c = —0.297952 . . . that was computed in Ref. [10] by enumerating a somewhat awkward 
double series obtained after a lengthy calculation. This constant c characterizing the leading finite size behavior is 
nonuniversal and is expected to depend on the details of the probability density f(x) of the noise. A natural question 
is: can one calculate this constant for arbitrary density function f(x)l In this paper we provide an exact formula for 
this constant c valid for arbitrary symmetric f(x). 

Our results are twofold. First we consider the class of density function f(x) that has a finite second moment, 
cr 2 = x 2 f(x)dx. Then a denotes the characteristic length of a single jump. Since, E[M n ] has the dimension of 
length, it is preferable to consider the dimensionless variable E[M n ]/cr. We show that for large n 



E(M n ) 
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The leading y/n behavior is universal (does not depend on the details of the density function f(x)) and easy to 
compute by appropriately rescaling the continuous time Brownian result. Our main new result is to obtain an exact 
expression for the nonuniversal constant 7. Our result is best expressed in terms of the characteristic function, 
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f{x) e lkx dx 



(5) 



For density functions with a finite second moment, i.e., when f(k) = 1 — a 2 k 2 /2 + 0(k ) as k — > 0, we show that 
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We prove in appendix- A that 7 < for arbitrary f(x), a fact not apriori obvious. Let us quote a few examples where 
the integral in Eq. (6) can be performed explicitly, 

(7) 
(8) 
(9) 
(10) 

where ((z) is the Riemann zeta function (analytically continued for z < 1). For the uniform distribution over [—1,1], 
our exact formula in Eq. (6) reproduces very simply the result obtained in Ref. [10]. In this case, using f(k) = 
and cr=l/ \/3 in Eq. (6) one gets 



m = 


\[8{x + 


1)+J(a:-1)] => 7 = -1/2 = 


-0.5 


m = 


I e -M 
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=> 7 = -1/V2 = -0.70710... 




m = 


a 2 


^ ^ 7= -(2V3- 1)/V6 = 


-1.00597 


m = 


1 


- x2/2 ° 2 => 7 = C(l/2)/v^ = 


-0.58259. 
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-0.516068... 



(11) 



where the integral was performed using Mathematica. Note that the constant c = —0.297952 ... in Eq. (3) is simply 
c = 70- with a — l/\/3 for the uniform distribution and 7 = —0.516068 . . . given in Eq. (11). 

We have also generalized our results to the case of Levy flights where the second moment diverges and one has 
f(k) — 1 — \ak\v + 0(k 2 ) as k — > with 1 < fj, < 2 [11]. Note that a has the dimension of length. The probability 
density f(x) of the step lengths has an algebraic tail for large |ar|, f(x) ~ In this case we show that for large 

n, the dimensionless expected maximum behaves as 



^1 = ^(1-1^^+7 + 0^-). 
a it \ [i J 



(12) 
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The leading term is again shown to be universal. We show that the leading finite size correction term is again a 
constant given by 



7 = 





[i -/(!)! 







(13) 



For example, for the case f(k) = exp[— \ak\^] with 1 < \l < 2, we obtain 
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The evaluation of the integral is presented in appendix-B. 

An interesting fact about the constant 7 is that even though it characterises the leading finite size correction to the 
expected maximum, it actually shows up even in the thermodynamic limit n —* 00 provided one looks at a behavior 
of a suitably defined quantity as follows. Let \x n \ denote the absolute value of the position of the walker after n steps. 
The distribution of x n , for arbitrary density f(x), can be computed relatively easily (see section-IV) and hence one 
can calculate i£[|a; n |]. We focus here on the case when the variance a 2 = x 2 f(x)dx as well as the fourth moment 
yU4 = x 4 f(x)dx of the jump distribution are finite. In that case, one can show that for large n, 
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3 ) 1+0(^-3/2). 

/ ^/n 



(15) 



Thus the leading term of i?[|a; n |]/cr for large n is exactly the same as that of the expected maximum E[M n ]/a in Eq. 
(4). However, unlike in the case of the maximum in Eq. (4), the leading finite size correction term in Eq. (15) is of 
0(n -1 / 2 ) and not a constant. Using Eqs. (4) and (15) one then gets 



E[M n ] - E[\x r , 



= 1 + 0{n- 1 ' 2 ). 



(16) 



Thus the difference between the expected positive maximum up to n steps and the absolute value of the expected 
final position of the walker after n steps, in units of a, tends to a negative constant 7 in the thermodynamic limit 
n — ► 00, a fact that is not apriori obvious. 

We end this section by mentioning a simple physical application of the results above. Let us consider the simplest 
model of a polymer chain namely the Rouse model [5] where the monomers are connected by harmonic springs. A 
configuration of the chain consisting of n monomers is specified by the position vectors {fi} of the monomers with 
i = 0, 1, 2 . . . , n. We assume that one end of the chain is grafted at the origin, rb = while the other end is free. We 
assume that the chain is at thermal equilibrium so that the probability of any given chain configuration is given by 
its Boltzmann weight, 



p [{n}} = ^- exp 



a n 
i=l 



(17) 



where Z n is the partition function, (3 — l/fc^T is the inverse temperature and the spring constant k characterises the 
harmonic coupling between neighbouring monomers. Let us now look at the components of the position vectors along 
any particular direction, say {x i \ 1 which can also be thought of as a one dimensional Rouse chain. The equilibrium 
weight in Eq. (17) indicates that the difference in position between the i-th and (i — l)-th monomer can be represented 
by a noise 



Xi-l 



(18) 



where &'s are independent and Gaussian distributed , /(£) = er^ l 2 ^ j o^FLk where a 2 — 1/(3k. Thus M n in Eq. (2) 
refers to the maximum displacement of the polymer chain along x direction and x n denotes the x coordinate of the 
end point of the chain (see Fig. 1). 
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FIG. 1. A typical configuration of the Rouse chain in 2-dimensions. M n denotes the positive maximum along the x direction 
and x„ denotes the x co-ordinate of the position of the end point. 

Thus, using the result in Eq. (10) in Eq. (16) we find that in the limit of a very long chain (n — > oo) at thermal 
equilibrium, the difference between the expected maximum displacement of the chain along a particular direction (say 
the x direction) and the absolute end to end displacement along the same direction tends to a nontrivial constant 



Note that in the context of the Rouse chain, the expectation E means a thermal equilibrium average over the Bolzmann 
weight in Eq. (17). The facts that (a) the difference approaches a constant and (b) that too, a negative constant, are 
not apriori obvious for the Rouse chain. 

The rest of the paper is organized as follows. In Section II, we set up the basic integral equation for the distribution 
of the maximum and provide an exact solution in Section II-A for the special case when the jump density is exponential. 
Section III deals with the general jump distribution where we present the Pollaczek-Spitzer formula and extract the 
finite size correction term exactly from an asymptotic expansion of this formula. These results are generalized to Levy 
processses in Section III-B. In Section- IV we calculate the expected value of the absolute position of the end point 
and demonstrate how the constant 7 shows up in the thermodynamic limit. We conclude with a summary and open 
problems. Explicit derivations of some of the formulae and integrations are relegated to the two appendices. 



In this section we set up an integral equation satisfied by the distribution of the maximum of a random walk for 
arbitrary symmetric jump distribution. We consider a random walk starting at x — at step n = and evolving 
via Eq. (1) where the noise £„'s are i.i.d. variables drawn from the common symmetric distribution Prob(£„ < x) = 
I- f(v) dy- We would like to compute the distribution of the maximum M n = max(0, X\,X2, ■ ■ ■ , x n ) up to n steps, 
i.e., Prob[M„ < y\. To derive this, we first define Q n {x,y) as the probability that, starting at x = x, the maximum 
of the walk up to n steps is less than or equal to y. Evidently, Prob[M„ < y] = Q n (0,y). Consider the first step 
where the particle jumps by an amount x\ — x which occurs with a probability density f(x\ — x) (see Fig. 2). 




= -0.58259.... 



(19) 



II. AN INTEGRAL EQUATION FOR THE DISTRIBUTION OF THE MAXIMUM 



y 



X 








n 



FIG. 2. A random walker, starting at x at n = 0, makes a flight to xi < y at step n = 1. 
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It follows, using the Markov property of the walk, that Q n (x,y) satisfies the following recursion relation, 

Q n (x,y) = ( Qn-\{xi,y) f(xi - x)dx 1 (20) 



with the initial condition Qo(x,y) — Q(y — x) where 6(z) is the Heaviside step function. Due to the translational 
invariance, it is also clear that Q n (x, y) depends only on the difference z — y — x, i.e., Q n (x, y) = q n (z = y — x) where 
z > 0. Making the change of variable, z' = y — x\ and using the translation invariance, the recursion in Eq. (20) 
becomes simpler, 

/•OO 

q n (z) = / qn-^fiz-z^dz', (21) 



valid for all z > and starting with qo(z) — 6(z). Thus, if one finds the solution q n (z) of Eq. (21), then the distribution 
of the maximum is just Prob[M„ < y] = Q n (0,y) = q n (y)- The density of the maximum is q' n (y) — dq n /dy. Hence, 
the expected maximum is E[M n ] = J °° q' n (y)ydy, the quantity we are after. 

The generating function, q(z, s) — Yl^Li Qn(z)s n then satisfies an integral equation 

/•oo />oo 

q(z,s) = s / q(z',s)f(z- z')dz' + s / f(z-z')dz', (22) 
Jo Jo 

valid for all z > 0. This integral equation is an inhomogcneous Wiencr-Hopf equation [12] and in general, for arbitrary 
kernel f(z — z'), it is very hard to solve this integral equation. The main source of difficulty is the fact that the limits 
of the integral on the right hand side of Eq. (22) are and oo, as opposed to say — oo and oo [12,2]. However, when 
the kernel f(z — z') is a normalized probability density function we may use the Pollaczck-Spitzer formula [22,21], to 
which we will come back to in section III. But, before that, it is instructive to solve Eq. (22) explicitly for special 
cases, whenever possible. In Section II-A, we solve Eq. (22) explicitly for the exponential density function. 



A. Exponential density function: An exactly solvable case 

For the exponential density functon, f(z) = ^e - ' 2 ', one can obtain an exact solution of Eq. (22). We first assume 
that the integral Eq. (22) is valid for all — oo < z < oo, though we are interested only in the solution for z > 0. Next, 
we note the identity, f"(z) = f(z) — S(z), where f"(z) — d 2 f/dz 2 . Diffcrcnting twice Eq. (22) and using the above 
identity one readily converts the integral equation into the following differential equation, 

^ = [l-ae(z)]q-a9(z). (23) 
For z > 0, the general solution is readily obtained, 

q(z, s) = + A(s)e-^~ sz + B(s)e^~ sz , (24) 

where A(s) and B(s) are two arbitrary constants (independent of z). Now, in the limit z — > oo, q n (z) — > 1 for all n, 
since the probability that the particle, starting at z — > oo, will not cross up to any finite step n is 1. Thus, one 
expects that as z — > oo, q(z, s) — > s/(l — s). Using this boundary condition in Eq. (24), one gets B(s) = 0. The other 
constant A(s) will be fixed by the matching conditions at z = 0. 
Now, for z < 0, the solution of Eq. (23) is given by 

q{z,s) = C{s)e z , (25) 

where we have used the boundary condition, q(z — > — oo, s) — » 0. Now, we are ready to match the solution in Eq. 
(24) for z > with that in Eq. (25) for z < 0. The continuity of the solution at z = and also the continuity of the 
first derivative at z = fixes the two constants A(s) and C(s) uniquely. We get, 

A(s) = — ; C(s) = — . (26) 

1 — 5 VI — S 

Thus, for z > 0, the exact solution of Eq. (22) is given by 
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The probability density for the maximum then has the generating function, 



£^ Wfl » = ^}Oe-^'. (28) 



n=l 

Subsequently, the generating function for the expected maximum is given by 

1 



^£[M„]s" = ^>" / q' n (z)zdz 

n=l n=l ^° 



(l_«)3/2 l_ s - 



(29) 



Expanding the right hand side of Eq. (29) in powers of n, we get 



£7[M n ] - -1 + _ r(n+i) , (30) 

where T(x) is the Gamma function. Note that the result in Eq. (30) is valid for all n > 0. The variance of the step 
lengths for the exponential density is given by, a 2 — z 2 e~' z 'cLz/2 = 2. Hence, the dimensionless number E[M n ]/a 
is given by 

E[M n ] _ 1 ; [2 T(n + 3/2) 



V2 V 7T T(n+1) 

f-J ? + 0( „-v, ) , (31) 

where we have made the asymptotic expansion for large n in the second line. Thus, the result in Eq. (31) is of the 
general asymptotic form as in Eq. (4) with the nontrivial constant 



III. THE GENERAL CASE: ASYMPTOTIC EXPANSION OF POLLACZEK-SPITZER FORMULA 

In this section, we derive the exact asymptotic behavior in Eq. (4) for the expected maximum for an arbitrary, 
symmetric jump density function f(z). As mentioned in the previous section, the solution of the integral equation 
Eq. (22) is hard to obtain analytically for an arbitrary kernel f(z — z'). However, when f(z) is a probability density 
function, Pollaczek derived a general formula giving the Laplace transform of the probability density of ordered partial 
sums of random independent variables [22]. In the special case of the distribution of the maximum, this formula was 
rederived by Spitzer [21] by a combinatorial approach. In principle, this solves the problem. However, extracting 
the precise asymptotic behavior of the first moment of the maximum is still nontrivial and this is what is precisely 
achieved in this section. Consider the following Laplace transform, 

E[e-" M «]= r e-"W n (z)dz, (33) 
Jo 

where q' n (z) = dq n /dz is the probability density of the maximum and q n {z) satisfies the recursion relation in Eq. (21). 
The Pollaczek-Spitzer formula for the generating function of the above Laplace transform [13] reads 



f^s n E[e-" M ^ =_L=0( a ,p) ; where 0( S , p) = exp 



In (l - sf{kj) 
p 2 + V dk 



(34) 



where < s < 1 and f(z) is an arbitrary symmetric normalized density function. In Eq. (34), f(k) — f{z)e lkz dz 
is the Fourier transform of f(z). 
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The generating function for the expected maximum can then be obtained by differentiation, 



To determine the asymptotic behavior of E[M n ] for large n, we need to know the behavior of h(s) near its principal 
singularity s = 1. It then follows from Eq. (35) that we need to know the precise behavior of the function </>(s, p) near 
s = 1 and p = 0. Below we analyse these asymptotic behaviors separately for two cases : (i) For density functions 
with a finite second moment, so that /(fc) — > 1 — a 2 k 2 /2 + 0(fc 4 ) as k — > where ct 2 is the variance of the jump 
lengths and (ii) for Levy flights where the jump lengths are power law distributed so that /(fc) — ► 1 — \ak\^ + 0(k 2 ) 
as k — > where a is a microscopic length and 1 < p, < 2. 



A. Jump lengths with a finite variance 



In this case, f(k) = 1 — a 2 k 2 /2 + 0(fc 4 ) as fc — > 0. To analyse 0(s, fc) near s = 1 and p = 0, it is first necessary to 
extract the most singular part of 4>(s, k) near s = 1 and p = 0. To do this, we first rewrite 



In (l - s/(fc)) = In M - s f 1 - \ 2 k 2 



In 



1 - S /(fc) 



1 - S (1 - ±CT 2 fc 2 ) 



(36) 



We next substitute Eq. (36) in the expression for (f>(s, p) in Eq. (34) and subsequently perform the first integral using 
the following identity [14] 



Jo 



ln(l 



s + sa 



2 k 2 /2) 



p 2 + k 2 



dk = ^ In (Vl - s + ap^fTJ 7 !^ . 



This gives 



<t>{s,p) 



exp 



dk 



In 



1 - sf(k) 



Vl — s + crp^s/2 
The extraction of the most singular part near s = 1 and p — gives 

1 



Wo P 2 + fc 2 \ l-s + scr 2 fc 2 /2 



(37) 



(38) 



cj)(s,p) 



[vT^7 + ap/V2\ 



exp 



^ Jo 



dk / 1 - /(fc) 
A; 2 I a 2 fc 2 /2 



where w in Eq. (39) means, in the strict mathematical sense, the following identity, 



lim — In 

s->l,p->0 p 



(VT - s + apv/^2) p)] = -i J 



'■ » rffc ln ( l - /(fc) 



fc 2 



(7 2 fc 2 /2 



Taking the derivative with respect to p and putting p = one gets from Eqs. (39) and (35), near s = 1 



h(s) 



V2(l-s) 3 / 2 7r(l-s)/ fc 2 



dk, 
ttt m 



2_ ^ 1 - /(fc) 



a 



fc 2 



(39) 



(40) 



(41) 



Noting that the singular behavior h(s) ~ (1 — s)~P (j3 > 0) of the generating function near s = 1 translates into 
the estimate E[M n ] ss n 13 ^ 1 /T((3) for large n, we get from Eq. (41) the following exact asymptotic behaviors of the 
expected maximum, 



2n 1 r dk. 



_2_ / 1 - /(fc) 
a 2 I fc 2 



+ 0{n 



(42) 



Dividing by ct yields our main result announced in Eq. (4) with an exact expression for 7 given in Eq. (6). As a 
check, we find that for the exponential density, f(z) = e~l z l/2, i.e., with /(fc) = 1/(1 + fc 2 ), the formula in Eq. (6) 
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yields 7 = —l/y/2, in agreement with the direct solution in Eq. (32). A few other cases where the integral in Eq. (6) 
can be performed explicitly are listed in Eqs. (7)-(10). For the case of uniform density, i.e., f(z) = 1/2 for — 1 < z < 1 
(and otherwise), we thus obtain an exact closed form expression for 7 whose numerical value is in agreement with 
the result obtained by Coffmann ct. al. by a different method [10]. Our result is evidently more general and holds for 
arbitrary symmetric jump density function. In appendix-A, we will prove that 7 given by Eq. (6) is always negative 
for arbitrary f(z). 



B. Levy distributed jump lengths 



We now consider the case when the second moment of the jump distribution diverges. In particular, we consider 
Levy jumps such that the Fourier transform of the density function behaves as f(k) = 1 — |aft|" + 0(k 2 ) for small k 
with 1 < p, < 2. This indicates that in real space the steps lengths have an algebraic tail, f(x) ~ |ir|~ 1_At for large \x\. 
Note that the Spitzer's formula in Eq. (34) is still valid for such processes. However, as we will see, the asymptotic 
behavior of E[M n ] is quite different from the case where a 2 is finite. 

We proceed as in the previous subsection by extracting the most singular behavior of 4>(s, p) near s = 1 and p = 0. 
We first rewrite, for k > 0, 



m (1 - ./(*)) = m (i - s (i - (aft)")) + m 



:(1 - {akY) 

We next substitute Eq. (43) in the expression for <j){s,p) in Eq. (34). This gives 

<j){s,p) = exp [-h(s,p) - h{s,p)] , 

where 

dk 



(43) 



(44) 



f 

h{s,p)=P- \ 

7f Jo 



p 2 +k 2 



In (1 - s + s(ak)») 



and 



-ln(l-,s) + 



>(«,p) = £ f 
* Jo 



i r ~ 

Wo 1 



du 



+ u 2 



In 



1 - s 



(apuY 



dk 



p 2 + k 2 



In 



1 - sf(k) 



1 - s + s(ak)t* 



(45) 



(46) 



Taking the derivative with respect to p and keeping only the leading singular terms near s = 1 and p = 0, we get 



d<t>{s,p ) 

dp 



dk 



7Ta/1 - s J a k 2 



In 



1 - sf(k) 
{aky 



pa^p^ 1 



i» du 



nVT^l Jo (1 +u 2 )(l - s + {apuYY 



(47) 



The second integral on the right hand side can further be simplified by first making a change of variable pu = y and 
subsequently taking the limit p — > 0. The resulting integral can be performed in closed form. Putting everything 
together, we find the leading singular behavior of h(s) near s = 1 and p = 0, 



h(s) = y 8 »E[M n ] = ^(//^i-y^) + _j_ r 

V ' <L^ 1 J ^1- e U+VM tt(1 - S) J 



n=0 



dk, 1 



sf(k) 



{akY 



+ 



1 



(1 - s y/n 



Tr(l-sY +1/ v 

where B(x,y) is the standard Beta function. Subsequently, one obtains the following large n behavior of E[M n ], 



(48) 



E[M n ] 



k 2 



{akY 



+ 0{n 1 ^- 1 ). 



(49) 



Note that for p = 2 and a — Eq. (49) reduces to Eq. (42), as it should. The dimcnsionless expected maximum 

is obtained by dividing Eq. (49) by the microscopic length a and one arrives at the result in Eq. (12) with the 
constant 7 given by the exact formula in Eq. (13). 
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IV. APPEARANCE OF THE CONSTANT 7 IN THE THERMODYNAMIC LIMIT 



In the previous section, we have demonstrated how the constant 7 appears as the leading correction term to the 
asymptotic ^fn behavior of the expected maximum E[M n ] for large n. In this section we show how 7 appears as the 
leading term in the large n limit if one considers the difference (E[M n ] — S[|a;„|]) /a, where £7[|a; n |] is the expected 
absolute end to end distance of the walker after n steps. 

The calculation of _E[|a; n |] is relatively straightforward compared to that of E[M n \. We start with the random walk 
in Eq. (1) where the jump density /(£) is symmetric and has a finite second and fourth moment, a 2 = J_ x 2 f(x)dx 
and j«4 = x 4 f{x)dx. Let P n (x) be the probability density for the particle to be between x and x + dx after n 
steps, starting from x = at n = 0. Using the Markov property of the walk, it is easy to see that P n (x) satisfies the 
recursion relation 



P n (x) 



f 

J —i 



P n (y)f(x-y)dy, 



(50) 



starting from the initial condition, Pq(x) = 5(x). Note that unlike the recursion relation in Eq. (21), the limits of 
the integral in Eq. (50) are respectively —00 and +00 and hence Eq. (50) can be easily solved by taking the Fourier 



transform, P n (k) = J? P n {x)e lkx . One gets P n (k) = f(k) 



from the inverse Fourier transform, 



Hence the solution, valid for all n > 0, is obtained 



P n (x) = - 



f(k) 



i k x 



dk. 



(51) 



The large n behavior of P n (x) can then be obtained by a straightforward scaling analysis of the integral in Eq. (51). 
In the large n limit, the most important contributions to the integral will come from the neighbourhood of k = 
where one can expand f(k) = 1 — a 2 k 2 /2 + ^ 4 fc 4 /24 + 0(fc 6 ). This gives, near k = 0, 



In 



f(k) 



' 2 k + 24 



3a 4 



k 4 + 0(k 6 ) 



Substituting in Eq. (51) and rescaling k^/n = q and y = x/y/n one gets 

Za 4 )q 4 



Pn{x) 



J_ r°° dq 

Vn J -co 2?r 



exp 



1 

sfn 



dq 
2n C 



24n 

(/M - 



3a 4 )q 4 



24n 



+ o( q y n 2 ) 

+ 0(q 6 /n 2 ) 



(52) 



(53) 



where we have expanded the exponential for large n holding y = x/^/n fixed. Using this probability density in Eq. 
(53), it is then easy to calculate the expectation 2£[|x n |] and one gets asn^oo, 



E[\x n \ 



2n 1 

12v/2tt 



3a 4 



-3/2 



), 



(54) 



in agreement with the result obtained by Petrov [23] using a somewhat different method. 

Comparing Eq. (54) with Eq. (42), one sees that the leading yfn term in both _E[|a; n |] and E[M n ] have the same 
coefficient. However, while the next subleading term in E[M n ] is a constant, the one in i?[|a;„|] decays as rT 1 / 2 for 
large n. Taking the difference of Eq. (42) and Eq. (54) and dividing by a, we obtain Eq. (16). Thus, as n — > 00, the 
difference between the two quantities (scaled by a) approaches a dimensionless constant 7 given by Eq. (6). Hence it 
is possible to observe the constant 7 even in the thermodynamic limit provided one looks at the difference between 
two observables. 



V. CONCLUSION 



In this paper, we have studied analytically the finite size corrections to the asymptotic large n behavior of the 
expected maximum S[M n ] of a one dimensional random walk of n steps with arbitrary, symmetric jump distribution. 
While the leading ^/n behavior is universal and easy to understand by extrapolating the result of continuous time 
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Brownian motion, the leading finite size correction term turns out to be a nonuniversal constant 7 which is nontrivial. 
In this paper, we have presented an exact formula for this constant, valid for arbitrary symmetric jump distribution. 
We have also generalized our results to the case of Levy processes. 

We have also demonstrated how this constant appears even in the n — > 00 limit as the leading behavior of the 
difference E[M n ] — E[\x n \], where x n is the position of the walker after n steps. As a nice application, we considered 
a Rouse polymer chain consisting of n monomers at thermal equilibrium and showed that the difference between the 
expected positive excursion along any direction (say the x direction) and the expected absolute value of the end to 
end distance of the chain approaches a negative constant in the thermodynamic limit 

where E denotes the thermal average and a = l/^/JTH (/? being the inverse temperature and k being the spring 
constant of the chain). The result in Eq. (55) is nontrivial and somewhat counterintuitive. 

There are several possible extensions of this work. In this paper, we have only considered the finite size behavior of 
the expected maximum. It would be nice to extend this finite size study to higher moments of the maximum, or even 
to the full distribution of the maximum. The distribution of the maximum of a set of correlated random variables is a 
subject of current interest and several papers have recently studied the distribution of the maximum or related objects 
for the continuous time Brownian motion, in the context of fluctuating interfaces [17-19] and also in the context of 
convex polygons and queueing theory [20]. For discrete step random walks, the continuous time results (suitably 
rescaled) will provide the leading asymptotic behavior of the distribution for large n. It would be interesting to see 
if the method presented in this paper can be used to investigate the finite size effect in the maximum distribution in 
these systems. 

It would also be interesting to extend our results to higher dimensions. For example, is there a leading constant 
correction term to E[M n ] for large n, where M n is the maximal radial distance from the origin of a random walker 
of n steps in d-dimensions? If so, it would be interesting to calculate this constant. For this, one needs to develop a 
higher dimensional analogue of the Pollaczek-Spitzcr formula which would be interesting in its own right. 



APPENDIX A: PROOF THAT 7 < 



In this appendix we show that the constant 7 that appears in Eq. (6) or Eq. (13) is negative for arbitrary symmetric 
jump density function f(x). For convenience, we first rewrite 7 = C/a where 



C* 



1 f°° dk 
Wo 



In 



1 ~ f(k) 

<7 2 k 2 /2 



(Al) 



where f(k) = J__° f{x)e tkx dx. We assume that f(x) is symmetric with a finite second moment a 2 . Clearly, the 
argument inside the logarithm in Eq. (Al) is positive. To prove that C < 0, it is sufficient to prove that the argument 
[1 — f(k)]/a 2 k 2 /2 inside the logarithm is less than 1. For symmetric f(x) one can write f(k) = 2 f(x) cos(kx)dx. 
The next step is to write the identity 



(l-/>) 
a 2 k 2 /2 



k 2 a 2 



/(*) 



cos(kx) — 1 



k 2 x 2 



dx, 



(A2) 



which can be proved by carrying out the integration on the right hand side explicitly and using the definition a 2 = 
2 Jp 00 x 2 f(x)dx. To prove that the left hand side of Eq. (A2) is less than 1 we need just to prove that the second term 
on the right hand side is positive. Using the elementary inequality, cos(kx) — 1 + k 2 x 2 /2 > for all x and the fact 
that f(x) > 0, it follows immediately that indeed the second term is positive. Thus one has 



(l-/(fc) 
a 2 k 2 /2 



< 1, 



(A3) 



which then proves that C < and hence 7 < 0. A similar proof can be easily constructed to show that for Levy 
processes as well 7 as given in Eq. (13) satisfies the inequality 7 < 0. 
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APPENDIX B: PROOF OF AN IDENTITY 



In this appendix we prove the following identity valid for all 1 < /j, < 2, 



1 f°° dk 
Wo 



In 



1 - e 



C(Vm) 



(2tt) 1 /m s in(7r/2^)' 



(Bl) 



where Q(z) — J2m=i m z 1S the Riemann zeta function which is usually convergent for z > 1. However, it is possible 
to analytically continue ((z) for z < 1 [15] and one defines 



C(z) = lim 



E 



1 - z 



We first make a change of variable fc M = x in the integral in Eq. (Bl). This gives 

1-e" 



— /° 

M7T Jo 



dxx- 1 - 1 ^ In 



Next we make one integration by parts and use the fact that fi > 1. This yields 



k Jo 



1 



1 



e x — 1 x 



(B2) 



(B3) 



(B4) 



Note that each integral on the right hand side of Eq. (B4) is separately divergent near x = 0, though their difference 
is convergent. To make progress, we introduce a small cut-off e and define 



K Jo 



1 



1 + e x + e 



- 1 ^dx = h{e)-I 2 {e). 



(B5) 



Eventually we are intersted in I = 1(0). The reason behind introducing this additional cut-off is so that the two 
integrals will be separately convergent for finite e and then after performing the two separate integrals, we will 
eventually take the e — > limit. 
The first integral can be written as 



1 r°° x- x i»di 

Wo e*-l + 



+ e 7r 



= _r l__ $ ( 1 - e,l ,1 



(B6) 



where $>(z,s,v) — J2^=o( v + n ) Sz " i s the Lerch function [14]. The second integral is elementary and can also be 
performed exactly to give 



h{e) = - 



1 f^x^/^dx e- 1 ^ 



x + e 



1 1 

Bl 1--,- 



where B(x, y) = T(x)T(y)/T(x + y) is the standard Beta function. Putting these together, we then have 



He 



r l- 



$ [ 1 - e, 1 , 1 



(B7) 



(B8) 



Now, the tricky part is to take the e — > limit in Eq. (B8). To do this, we will make use of the following asymptotic 
behavior of <f>(z, s, v) as z — > 1 [16] 



$(z, s, «) = T(l - s)[- lnz] 8 - 1 ^-" + C(s, a) 



(B9) 



where £(s, a) = X)^Lo( n + a ) ' '■ Putting z = 1 — e m Eq. (B9) and expanding for e, one then finds that as e — > 0, 
the two leading order terms are given by 



#(i-«,i-i.i)-r(i)^. + f( 1 _i, 1 ) +0(<! .-./» ) . 

Substituting this behavior in Eq. (B8) we get 



(BIO) 
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/./(«^o)-ir(i-i)c(i-i.i)-ir(i-i)c(i-i 



(Bll) 



At / \ A* / 77 V AU 

where we have used the fact that £(z, 1) = ((z). One can further rewrite Eq. (Bll) by using the identity [14], 

C(z) = 2{2tt) z - 1 sin(7rz/2)r(l - z)C(l - z). (B12) 

Using this identity in Eq. (Bll) one readily arrives at the final result in Eq. (Bl). 

In particular, note that the result in Eq. (Bl) is valid even for [i = 2. In that case one gets, I(fi = 2) = C(l/2)/ v / 7r. 
Note that when the jump lengths are Gaussian distributed as in Eq. (10), one has f(k) — exp[— k 2 a 2 /2]. Hence, from 
Eq. (6) 



7 = 



7T 



V2 



dk 



In 



k 2 



IJjm = 2) _ C(l/2) 
V2 



2tt 



= -0.58259. . 



(B13) 



thus proving the result in Eq. (10). The numerical value of ((1/2) 
using Mathematica and was used in Eq. (B13). 



-1.46035 can be obtained to arbitrary precision 



Note added in proof: We thank R.M. Ziff for pointing out that the constant c = ■ycr = —0.297952.. for the uniform 
case also appeared in an apparently unrelated three dimensional trapping problem first studied in [Ziff R. M., flux to 
a trap, 1991 J. Stat. Phys. 65 1217] where it was evaluated by numerically iterating a set of recurrence relations. 
Our Eq. (11) provides an exact expression of this constant. 



[1] S. Redner, A Guide to First-Passage Processes (Cambridge University Press, Cambridge, 2001). 

[2] W. Feller, An Introduction to Probability Theory and Its Applications (Wiley, New York, 1968), vol-II, p-381. 

[3] B. Hughes, Random Walks and Random Environments (Clarendon Press, Oxford, 1996). 

[4] D.E. Koshland, Bacterial Chemotaxis as a Model Behavioral System (Raven, New York, 1980). 

[5] P.E. Rouse, A theory of the linear viscoelastic properties of dilute solutions of coiling polymers, J. Chem. Phys. 21: 
1272-1280 (1953). 

[6] E.W. Montroll and M.F. Shlesinger, The wonderful world of random walks, in Nonequilibmum Phenomena II: From 
Stochastics to Hydrodynamics, Studies in Statistical Mechanics, vol-11, eds J.L. Lebowitz and E.W. Montroll (North 
Holland, Amsterdam, 1984). 

[7] A. Blumen, J. Klafter and G. Zumofen, in Optical Spectroscopy of Glasses ed. I. Zschokke (Reidel, Dordrecht, 1986), p-199. 
[8] J. P. Bouchaud and A. Georges, Anomalous diffusion in disordered media: statistical mechanics, models and physical 

applications, Phys. Rep. 195, 127-293 (1990). 
[9] E.G. Coffman and P.W. Shor, Packing in two dimensions: asymptotic average-case analysis of algorithms, Algorithmica, 

9: 253-277 (1993). 

[10] E.G. Coffman, P. Flajolet, L. Flato, and M. Hofri, The maximum of random walk and its application to rectangle packing, 

Probability in Engineering and Informational Sciences, 12: 373-386 (1998). 
[11] J. Bertoin, Levy Processes (Cambridge University Press, Cambridge, 1998). 

[12] P.M. Morse and H. Feshbach, Methods of Theoretical Physics (McGraw-Hill, New York, 1953), vol-I, p-978. 
[13] F. Spitzer, The Wiener-Hopf Equation whose kernel is a probability density, Duke Math. J. 24: 327-343 (1957). 
[14] I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, series and Products (Academic, San Diego, 1994). 
[15] M. Abramowitz and LA. Stegun, Handbook of Mathematical Functions (Dover, New York, 1973). 

[16] W. Magnus, F. Oberhettinger and R.P. Soni, Formulas and Theorems for the Special Functions of Mathematical Physics, 

(Springer- Verlag, Berlin-Heidelberg, 1966), p-33. 
[17] S. Raychaudhuri, M. Cranston, C. Przybyla, and Y. Shapir, Maximal height scaling of kinetically growing surfaces, Phys. 

Rev. Lett. 87:136101 (2001). 

[18] S.N. Majumdar and A. Comtet, Exact maximal height distribution of fluctuating interfaces, Phys. Rev. Lett. 92: 225501 
(2004); see also cond-mat/0409566 (to appear in J. Stat. Phys. (2005)) for a longer version. 

[19] G. Gyorgyi, P.C.W. Holdsworth, B. Portelli, and Z. Racz, et. al., Statistics of extreme intensities for Gaussian interfaces, 
Phys. Rev. E 68: 056116 (2003). 

[20] M.J. Kearney and S.N. Majumdar, On the area under a continuous time Brownain motion till its first-passage time, J. 
Phys. A: Math. Gen. 38: 4097-4104 (2005). For a discussion of the discrete perspective of queues, see M.J. Kearney, On 



12 



a random area variable arising in discrete-time queues and compact directed percolation, J. Phys. A: Math. Gen. 37: 
8421-8431 (2004). 

[21] F. Spitzer, A combinatorial lemma and its application to probability theory, Trans. Am. Math. Soc. 82: 323 (1956). 
[22] F. Pollaczek, Fonctions caracteristiques de certaines repartitions definies au moyen de la notion d'ordre, Comptes rendus 
234: 2334 (1952). 

[23] V.V. Petrov, Sums of Independent Random Variables (Springer, Berlin, 1977), chapter 8. 



13 



